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We calculate the strength of symmetry restoration effects in highly 
non-equilibrium states which can arise, for example, during preheating after 
inflation. We show that in certain parameter range the one-loop results are 
unstable, requiring summation of multiloop diagrams. We solve this problem 
for the O(N) model in the large iV-limit and show that the symmetry restora- 
tion may be less effective than what predicted by the one-loop estimate. 
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It has been recently realized [T|] that the process of reheating after the inflationary 
era || may consist of different epochs. 

In the first stage, the effects of stimulated decays and annihilation may lead to an 
extremely effective explosive particle production. The energy is released in the form of in- 
flaton decay products, whose occupation numbers are extremely large, and have energies 
much smaller than the temperature that would have been obtained by an instantaneous 
conversion of the inflaton energy density into equilibrated radiation. 

Since it requires several scattering times for the low-energy decay products to form 
a thermal distribution, it is rather reasonable to consider the period in which most of 
the energy density of the Universe was in the form of non-thermal quanta produced by 
inflaton decay as a separate cosmological era. Following we refer to the final stages 
of parametric resonance as the preheating epoch to distinguish it from the subsequent 
stages of particle decay and semiclassical thermalization || when the inflaton decay can 
be still efficient, but amplitudes of fluctuating fields became gradually smaller because 
of Universe expansion. 

Several aspects of the theory of explosive reheating have been studied in the case of 
slow-roll inflation @] and first-order inflation ||. In particular, the phenomenon of sym- 
metry restoration during the preheating era has been investigated recently by Kofman, 
Linde, and Starobinski || and by Tkachev j7j in the framework of typical chaotic infla- 
tionary models. It was shown that nonthermal symmetry restoration processes during 
the nonequilibrium stage of preheating may be very efficient with important implications 
for models containing relatively low scales like the invisible axion model or supersym- 
metry with Polonyi fields, and may be for Grand Unified Theories (GUT). For example, 
if a GUT symmetry is restored during the preheating epoch, the subsequent symmetry 
breaking phase transition will reintroduce the problems of monopoles M or domain walls 
||. In any case, if one askes whether particular high-energy symmetries are restored after 
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inflation or not, one has to study the non-equilibrium preheating epoch. 
We consider a model with the following potential 

W,r ? ) = A0 4 + ^(r ? 2 -r ?o 2 ) 2 + |0V- (1) 

The field <fi corresponds to the inflaton. While inflaton self-coupling is restricted to be 
A 10 -13 , there is no such severe restrictions on other couplings. Among products of 
the inflaton decay there can be different species, but for simplicity we assume here that 
quanta of the field r\ constitute the major channel of its decay, which means in particular 
that g ^> A. We assume also that 77 transforms as a vector under the action of O(N), 
i.e. rf = r] a 7] a with a — 1, ■ ■ -, N. 

In previous papers [0,01,0 the effects of the coupling a was neglected, which corre- 
sponds, in fact, to the one loop approximation. The importance of multiloop correction 
was noticed in Ref. J7|. We shall show here that they change the results dramatically 
when a ^> g and the one-loop approximation actually breaks down in this region. 

We do not consider an implications of non-zero a on parametric resonance itself, but 
consider the effects of symmetry restoration only. We assume that at the end of the 
resonance stage the energy density which builds up in the fluctuations of field 77, and 
which we denote as p v , is some fraction of inflaton energy density. Simple assumption 
would be that this fraction equals to one half |3J . Numerical integration of Ref. for 
the case of simplest A0 4 model confirms that the resonance ends when half of the energy 
is in fluctuations, however it might not be true for g ^> A. Morover, because of the 
expansion of the Universe during preheating, p v is many orders of magnitude smaller 
than the initial inflaton energy density. Since this is model dependent, at present we 
keep free parameter. 

We can parametrize the distribution function of the created quanta at preheating as 
/ = f(k/k*) where k is particle momenta. What is important here is the smallness of 
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typical particle momenta, k*, with respect to the temperature that one would get in the 
case of istantaneous reheating neglecting the expansion of the Universe. In the extreme 
case when all inflaton couplings are not very different, fc* is of order of inflaton mass, 
fc* ~ 10 13 GeV. For defmitness, let us assume that the distribution function is of the 
form 

f(k) = A8(\k\-\K\). (2) 

Parametrically, and that is what we are interested in here, Eq. (Q) will reproduce the 
correct unswer. Moreover, the distribution is really peaked at the end of the resonance 
stage in reasonably wide range of model parameters. 

The constant A (or its initial value in our problem) can be fixed by computing the 
energy density of 77 particles and setting it equal to p v . Notice that, since the order 
parameter changes in time, e.g. during the process of symmetry restoration, the energy of 
particles coupled to it changes too, but three-momenta keep constant in an homogeneous 
background. Moreover, in this case the effective potential simply coincides with the 
energy density, which simplifies the calculations. To answer the question whether the 
symmetry tends to be restored or not, it is sufficient to consider homogeneous background 
only. 

As opposed to large-angle scattering processes, forward-scattering processes do not 
alter the distribution function of the particles traversing a gas of quanta, but simply 
modify the dispersion relation. This remains true also in the case of a nonequilibrium 
system. Forward scattering is manifest, for example, as ensemble and scalar background 
corrections to the particle masses. Since the forward scattering rate is usually larger 
than the large-angle scattering rate responsible for establishing a thermal distribution, 
the nonequilibrium ensemble and scalar background corrections are present even before 
the initial distribution function relaxes to its thermal value. These considerations allow 
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us to impose k° = y |k| 2 + m 2 (rj) as the dispersion relation for the particles created by 
the parametric resonance, where fj is the order parameter, with fj = tjq in vacuum, and 
m 2 (fj)5 ab = a [(fj 2 - r)l)8 ab + 2fj a fj b ] /3!. 

Let us first consider the case N = 1. We can not use the imaginary-time formalism 
to determine the effective potential for the scalar field fj during the nonequilibrium pre- 
heating period since in the nonequilibrium case there is no relation between the density 
matrix of the system and the time evolution operator, which is of essential importance 
in the formalism. There is, however, the real-time formalism of Thermo Field Dynamics 
(TFD), which suites our purposes ||10|| . This approach leads to a 2 x 2 matrix structure 

for the free propagator (only the (ll)-component is physical) 

( D n (K) D 12 (K) \ / A(K) W f(k) 6(k ) + f(k) \ 

\D 21 (K) D 22 (K) ) \ A* (if) ) \ 0{-k o )+f(k) f(k) ) 

x 27:S[K 2 -m 2 rl (fj)], (3) 
with the usual vacuum Feynman propagator 

The one-loop effective potential Vi(fj) can be easily obtained solving the equation 

dVi ^ If d A K n /r , N 1 , r 

T i = o / 77^1 D U (K)^-(I + I f ), (5) 



dm 2 1 2 J (2vr) 4 i±v ' 2 
where Ti is the one-loop tadpole diagram which receives contribution only from the 

diagram in Fig. 1 (since the final external leg is fixed to be that of the physical field) 

and we have separated the usual zero-temperature contribution from the one given by 

the gas of //-quanta with distribution function /, Eq. (2). One then obtains 

V, (fj) = Vl 0) + V{ f) = Vf } + y/\K\ 2 + m*(fj), (6) 

where is the one-loop zero-temperature effective potential. With the distribution 
function Eq. @, the particle number density is 



and the particle-dependent part of the effective potential reads v}^ = 72^/|k*| 2 + m 2 (r])/2. 
We would like to mention for the future use that Eq. (^) is valid only in the one-loop 
approximation . Expanding around f] = 0, we obtain that the coefficient of f/ 2 -term 
in the potential (or the effective mass of the field) is changed by interaction with the 
medium by an amount cm/ |k*| ~ ap^/lk*! 2 ||, J7|. The presence of nonthermalized 
77-quanta may lead to symmetry restoration if p,j/|k*| 2 > rj 2 . 

Our main concern in this Letter is on what happens when multi-loop corrections 
are included in the computation of the nonequilibrium effective potential. We will show 
that in certain range of parameters the effective potential gets large contributions at the 
two-loop level and higher orders of perturbation theory, even larger than the one-loop 
contribution. Perturbation theory is then lost unless a proper resummation can be done. 
We give an example in which the resummation can be performed exactly. 

Let us first consider the theory described by the potential in Eq. (|l]). By applying 
the rules of TFD we obtain the two-loop tadpole diagrams as drawn in Fig. 2 (plus 
counter-terms arising at one-loop). Underlined diagrams are identically zero since they 
contain 5(K 2 - m 2 )5((K - P) 2 - m 2 )5(P 2 - m 2 ) = 0. 

We first consider the simplest case of diagram 2a 

J i I % I $ ^ D ^ K) ~ D ^ P) D ^ P) D ^ ■ (8) 

Observing that the tadpole does not have any imaginary part and that Re Dn(K) = 
KeD22{K), one can factorize the .fT-integral out and, using the mass- derivative formula 



TT|1 , obtain 

dV {2a) ia r d A K „ „ d 4 P 



I i^ ReD ^ K) i^ [A 2 (m/(A 2 (p)-(A* ( p)) 2 ) 



dm 2 4 J (2tt) 4 11V ' (2? 



" ! ' ( ' K ReD n (K) 



[^A 2 (P)- i ^L- 
J (2tt) 4 1 ] dm 2 



4 J (2vr) z 



where we have defined Jo — (2vr) -4 / d A P A 2 (P). Eq. (|9]) shows the absence of pinch 
sungularities, or A(P)A*(P) terms in the final result. The cancellation procedure works 
essentialy in the same way as in the equilibrium case and is totally independent of 
the distribution / |Tj|. After renormalization, the term proportional to If Jo gives a 
contribution ~ alflnm 2 which is suppressed by 0(a) with respect to the one- loop result. 
One should not claim victory to soon, though. Let us extract the /-dependent part only: 
its contribution goes like a If (dlf/dm^j. Expanding such an expression around fj — 0, 
we see that the contribution of diagrams of Fig. 2a is of order of (ap^/lk*) 4 ) with respect 
to the one-loop result. Since a may be ~ 1 and |k*| 4 /p^ <C 1 at the preheating era, we 
discover that the contribution from the two-loop tadpole diagrams is much larger than 
the one-loop result by several orders of magnitude and perturbation theory is lost!. 

Things get even worse when we consider p-loop diagrams. The tadpole diagrams 
shown in Fig. 3 (p > 2) contribute 



dm 2 2 V 2 / \ dm 2 ; 

All the terms which contain a vacuum contribution are ultraviolet divergent and by 
cancellation with the counter-terms we are left with only the /-dependent part which 
has the following behaviour (with respect to the one- loop result) 

(■2C-(»f" 

It is possible to add the tadpoles in a different way. Let us consider the p-loop 
contributions shown in Fig. 4, the so-called daisy diagrams. Using again the mass- 
derivative formula, we obtain the expression for a product of p > 2 successive propagators 
dV {4) _ 1 fa\P , x „ 1 / d x '" 



<K 2V27 + (/ » + / '»- (12) 

This expression has also to be renormalized. For p > 2 the ultraviolet singularities are 
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resummed in Jo and are cancelled by adding a series of counter-terms arising at one- loop. 
The result is 



It is clear that successive derivatives with respect to the squared mass give increasing 
inverse powers of the squared mass. The derivatives of Jo give a subleading contribution 
with respect to the /-dependent part. Extracting the /-dependent part and expanding 
such it around fj = 0, we discover that the contribution of the p-loop daisy diagrams is 
again of order of (a p r) /|k*| 4 ) p 1 with respect to the one-loop result. 

The discussion above shows that in order to obtain a more accurate information 
about the issue of nonthermal symmetry restoration one must study an infinite series of 
diagrams in perturbation theory. This is exactly analogous to what happens in a simple 
Xip 4 theory in equilibrium at finite temperature where the leading contributions to the 
effective potential in the infrared region come from the daisy and superdaisy multiloop 
graphs Hl3| . 

To deal with this problem we need a self-consistent loop expansion of the effective 
potential in terms of the full propagator. Such a technique was developed some time ago 
by Corwall, Jackiw and Tomboulis (CJT) in their effective action formalism for composite 



operators |L|J and may be also applied to nonequilibrium phenomena [15|]. One considers 
a generalization r[?7, G] of the usual effective action, which depends not only on fj(x), but 
also on G(x,y), a possible expectation value of the time-ordered product (T r)(x) r)(y)). 
The physical solutions satisfy the stationary requirements 
6T[fj,G] 



5rj(x) 



0. 

0. (14) 



5G(x,y) 

The conventional effective action r[^] is given by F[fj, G] at the solution Go(fj) of Eq. 



(|I4|). In this formalism it is possible to sum a large class of ordinary perturbation- 
series diagrams that contribute to the effective action T[fj], and the gap equation which 
determines the form of the full propagator is obtained by a variational technique. 

We now apply the CJT formalism in the limit of large N when the next-to-leading 
terms can be exactly summed. In each order we keep only the term dominant in N for 
large values of N. This allows us to resum exactly the series of the leading multiloop 
diagrams and to solve the gap equation for the full propagator without any approxima- 
tion. 

In order to obtain a series expansion of the effective action, one introduces the func- 
tional operator 

S 2 1 

D ab\v,x,y) = , (15) 

5r) a (x)5r) b (y) 



where I is the classical action. The required series obtained by CJT is then [[L4 

r[r/, G] = /(//) + ^ Tr In A) G^ + ^Tr [iT 1 G - l] + V 2 % G], (16) 

where Dq 1 = - (d^ + m 2 ) 5 ab 5\x, y), and T 2 [fj, G] is a sum of all two-particle- irreducible 
vacuum graphs in the theory with vertices defined by the classical action with shifted 
fields I[fj + rj\ and propagators set equal to G(x, y). 

Previous calculations show that among the multiloop graphs contributing to the 
effective potential in the 0(iV)-theory, only the daisy and superdaisy diagrams survive 
in the limit of large N |TJ|. This enables us to consider in I^^G] only the graph of 
0(a) given in Fig. 5. This is the Hartree-Fock approximation which is known to be 
exact in the many-body version of our large N limit. 

It is more convenient to concentrate on the effective masses rather than on the ef- 
fective potential. By stationarizing the effective action T[fj, G] with respect to G a t, we 
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obtain the gap equation 

G ab( x , V) = D dbi x , V) + T Kb G cc (x, x) + 2 G ab (x, x)] <5 4 (x, y). (17) 

o 

This equation is exact in the limit of large N and contains all the informations about 
the dominant iV-contributions to the full propagator. Indeed, the exact Schwinger-Dyson 
equation reduces to Eq. (|17D for large N. We Fourier-transform Eq. ( |P7D and take fj a = 0. 



The gap equation then reads [£L4 



M 2 = m 2 + — I AM 2 ) = m 2 + aN f , ^ . (18) 

We see that the one-loop results are stable when Ik*! 2 ^> M 2 , which translates to the 
condition p v <C |k + 1 4 /«, where now p v indicates the total energy of the noninteracting 
gas (the sum over all r] a ). In the opposite case the gap equation is approximately solved 
by M 3 ~ a N A |k*| 2 . Using p n ~ a N 2 A 2 |k*| 4 /M 2 we find 

M 2 ~ v /7i^. (19) 

In other words, the strength of symmetry restoration measured in terms of the effective 



temperature T eff /12 =< (rj — fj) > is given by T cff /12 « M /a ~ J p v /a. This result 
and Eq. (|T^) are intuitively understandable. Indeed, in this regime the contribution 
to the energy from self-coupling is important and p v is saturated by the self-interaction 
term in Eq. (1). This can be obtained in Hartree-Fock approximation directly applied 
to equations of motion for the 77-field (in the large iV-limit Hartree-Fock approximation 
becomes exact). 

To summarize. In the limit |k*| 4 3> ap v the one-loop results are stable and T 2 S ~ 
Prj/lk^ 2 . Note that our parameters, |k*| and p n are outcome of the stage of resonant 
decay and depend on the couplings g and A in Eq. (1). For example, we can expect 



that Ik*! 2 ~ y/g~% Mpj 0. In the opposite limit, |k*| 4 <C ap v , we have found T 2 t 



cfT ~ 
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^Jp v /a. This is smaller then the one-loop result by a factor |k*| 2 / yjotp^. This result 
remains valid when the inflaton decay products 77 and the order parameter in question 
correspond to different fields. The strength of the symmetry restoration in a highly 
non-equilibrium state cannot be traced by the one-loop result in the case of sufficiently 
strong interaction of inflaton decay products. This can have important consequences on 
the issue of symmetry restoration of various simmetries during the preheating era. 
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